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Abstract 

Takahashi & Soda (2010 Class. Quantum Grav. 27 175008) have recently con- 
sidered the effect (at lowest non-trivial order) of dynamical, quantized gravita- 
tional fluctuations on the spectrum of scalar Hawking radiation from a collapsing 
Schwarzschild black hole. However, due to an unfortunate choice of gauge, the 
dominant (even divergent) contribution to the coefficient of the spectrum cor- 
rection that they identify is a pure gauge artifact. I summarize the logic of their 
calculation, comment on the divergences encountered in its course and comment 
on how they could be eliminated, and thus the calculation be completed. 

PACS numbers: 04.70.Bw, 04.62. +v, 04.70. Dy, 05.40.-a 



1 Introduction 



In their recent work 10 , Takahashi & Soda have tackled an interesting and challenging 
question, that of assessing the influence of dynamical, quantized metric fluctuations 
on the spectrum of scalar Hawking radiation from a black hole. A natural sibling 
question that could be directly attacked with essentially the same techniques is that 
of back reaction of Hawking radiation on the quantum geometry of the black hole. 
An answer to either of these questions would provide us with valuable insight into the 
properties of quantum gravity, as seen through the prism of its reduction to effective 
quantum field theory on a curved background. 



The authors of 10 show that the spectrum of scalar Hawking radiation from a 
collapsing Schwarzschild black hole is corrected due to the presence of dynamical, 
quantized gravitational fluctuations interacting with the quantum scalar field at the 
lowest non-trivial (cubic) order. The choice of a collapsing spacetime instead of an 
eternal black hole is similar to the choice made in Hawking's original calculation J3j. 
The same result could be obtained by using the Unruh vacuum on an eternal black hole 



spacetime 11 . The corrected expected number of quanta (expected by an asymptotic 
observer at future null infinity) in a single field mode of static frequency uj and spherical 
harmonic index I, according to equations (71), (81) and (169) of 10 , is 



= + coth ( ™/ K) + • • • , (1) 



where Ci is a numerical constant, M is the Schwarzschild mass, k = 1/(4M) is the 
surface gravity, l v is the Planck length, e > is a small distance cutting off radial in- 
tegration before the horizon, and w cu t and L are parameters regulating the divergence 
of a 1-loop Feynman integral. Higher order corrections are presumed to be subleading 
in l p or one of the regulator parameters. 

I will summarize some of they key steps in the calculation leading to ([T]) and dis- 
cuss the divergences that appear along the way, as well as ways of resolving them. 
Section [2] reviews some relevant background material and establishes the notation. 
Section [3f discusses the divergence that prompted the introduction of the e- regulator, 
how it is related to the authors' choice of gauge and how a better choice of gauge 
eliminates this divergence and the need for this regulator. Section [4] discusses the di- 
vergence that prompted the introduction of the uj cut and L regulators and how these 
regulators could be removed using standard perturbative renormalization. Finally, 
section [5] emphasizes that, as a consequence of the results of section [3j the dominant 



contribution, identified in 10 , to the coefficient in front of the correction in equa- 
tion ([T]), is a pure gauge artifact. It also summarizes the non-trivial steps achieved 
in the calculation of Takahashi & Soda and how it could be completed to obtain a 
reliable, parameter-free estimate of the size of the correction to the Hawking spectrum. 



2 Background and notation 



The calculation in 10 starts out by quantizing free scalar and metric perturbation 
(graviton) fields on a collapsing Schwarzschild background. The Bogolubov coeffi- 
cients, which transform between asymptotic in- and out-modes, are estimated using a 
geometric optics approximation and give the standard Hawking spectrum for the free 
scalar field. Then, the metric perturbations are gauge fixed and reduced to two phys- 
ical scalar degrees of freedom, which are quantized as free fields. Finally, the explicit 
form of the interaction vertex between the scalar field and the metric perturbations is 
identified and used to compute the correction to the Hawking spectrum of the scalar 
field at lowest perturbative order. The relevant details of these steps are discussed in 
the rest of this comment. 

An important step in quantizing a free field is solving its classical equations of 
motion. Usually, this is accomplished by decomposing an arbitrary solution into a 
set of modes. On a Minkowski background, due to the translational symmetry, it is 
easiest to work with the set where each mode function is proportional to a plane wave 

4(x) = exp^x"), (2) 

where fc^ is the wave vector indexing the mode and x M are global inertial coordi- 
nates. We restrict ourselves only to the case of massless fields, hence fc 2 = 0. On a 
curved background, we can no longer make use of translation symmetry. However, the 



background used in 10 is presumed to be (at least outside the horizon) spherically 
symmetric, static in the asymptotic future, and approximately static in the asymptotic 
past. To make maximal use of the available symmetry, it is convenient to decompose 
the fields into a set of modes, where each mode function is proportional to a scalar 
or tensor spherical harmonic as well as 4>^^j{t,r) or 4>^j(t, r), where j indexes inde- 
pendent tensor polarizations (if any). These functions have the following asymptotic 
properties 

i?± (r) fart->±oo, (3) 



uutj e 
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where, up to normalization, the radial functions R^gj are uniquely specified by the 
geometry and the above listed conditions. The fact that the 4>~^j and mode 
functions do not coincide is ultimately responsible for the Hawking effect [2j Sec. 10.2]. 

While the decomposition of a scalar field into spherical harmonics is a standard 
exercise, the case of metric perturbations is more subtle, but is a topic with an ex- 
tensive literature, starting with the seminal work of Regge & Wheeler |7j. An up to 
date review can be found in [lj. Reference [5j is particularly useful as it presents the 
formalism of metric perturbations in a spherically symmetric spacetime in a way that 
is gauge invariant and covariant with respect to changes of coordinates in the (t, r)- 
plane. For reference, we establish the correspondence between the notations of [lO] 
and [5j. 

Due to spatial inversion symmetry, the perturbations naturally decompose into 
odd and even parity sectors. We shall only consider the even parity ones. Note that, 
below, asterisks denote components that can be deduced from symmetry. 

ffH H x v\ a \ 

Takahashi & Soda [To): h, w =\ * H/f w\ a (4) 

V * * r 2 Kj ab + B\ ab ) 

The components of the metric perturbation h^ v are given in static Schwarzschild 
coordinate^] /(r) = 1 — 1M.fr, ^ab is the standard metric on the unit 2-sphere, and 
T"'\ a ... denotes covariant differentiation of the tensor T with respect to "f a b- The 
components of the perturbation are parametrized by the scalars H, Hi, H, v, w, K 
and B. The components fH and H/f correspond to hu and h rr respectively. 

Martcl & Poisson (§]: p MV = ^ ( ab ^ r 2 K tm Q A BY emB + r ^G im Y im ) ^ 

tra AB / 

The components of the metric perturbation p p „ are given with respect to coordinates 
that respect spherical symmetry but are arbitrary on the (t, r)-plane, I and m are re- 
spectively the orbital and magnetic spherical harmonic indices and Qab is the standard 
metric on the unit 2-sphere. The components of the perturbation are parametrized 
by the scalars K im , G lm , and the (t, r)-plane tensors hlg and jf m . The vector and 
tensor spherical harmonics are defined from the scalar Y im as follows: Y^" 1 = D^Y 
and Y^§ = [DaD b + + l)FlAB]Y em , where Da is the covariant derivative with 
respect to SIab- Where no confusion is possible, the spherical harmonics indices may 
be omitted. 

The correspondence between the two notations should now be clear, given the 
equality h^dx^dx 11 — p^dx^dx 1 ' . Though, note that [lO] uses lowercase Latin in- 
dices, T a , for spherical tensors, while [5] uses lower case Latin indices for (t, r)-plane 
tensors and upper case Latin indices, Ta, for spherical tensors. Formulas given while 
referring to a given paper will use the notation of that paper. 



3 Singularity of 'convenient' gauge 

For canonical quantization of the metric perturbations, it is necessary to isolate the 
single physical degree of freedom in the even parity sector, known in the literature 

1 More precisely, as described previously, they are approximately static coordinates on the collaps- 
ing spacetime. 
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as the Zerilli or Zerilli-Moncrief function, the gauge freedom arising from linearized 
coordinate transformations must be completely fixed. Takahashi & Soda impose what 
they call a 'convenient' gauge, where v = 0, B = and K = 0. Further, they 
solve the constraints following from the equations of motion and express all remaining 
perturbation components in terms of the scalar Zerilli function ip z . I will show that 
gauge transformations required to enforce the 'convenient' gauge are singular at the 
horizon. As a consequence, some remaining non-zero components will be singular at 
the horizon as well. This singularity becomes obvious when these components are 
explicitly expressed in terms of the Zerilli function in a coordinate system regular at 
the horizon. 

Let a field be called regular at a point x if it is continuous and smooth in some 
neighborhood of that point; otherwise it is called singular at x. Regularity at x im- 
plies that the components of the tensor field must be continuous, smooth functions 
in a neighborhood of x, when expressed in any coordinate system that covers x. The 
converse implication is also true. It is a sensible mathematical and physical require- 
ment that metric perturbations be restricted to everywhere regular tensor fields. Note, 
however, that regularity in static Schwarzschild coordinates is insufficient to establish 
global regularity. Recall that Schwarzschild coordinates consist of two coordinate 
charts (exterior, r > 2M, and interior, < r < 2M), neither of which covers any 
point on the horizon. To check regularity at the horizon, it is necessary and suffi- 
cient to check continuity and smoothness in any coordinate system that does cover 
the horizon. We shall use the advanced Eddington-Finkelstein (EF) coordinates [4jj6]; 
they are regular on the future horizon, which suffices for our purposes. These (v, r) 
coordinates are related to the static Schwarzschild (f, r) coordinates as follows: 



dr 

— = r + 2Mln 



r 

2M ~ 1 



(6) 

v = t + r* dv = dt + /"Mr d t = d v 

r = r dr = dr d r — + d r 



3.1 'Convenient' gauge in EF coordinates 

The gauge chosen by the authors of [lO] for even-parity perturbations is termed the 
'convenient' gauge: K = B = v = 0. In the notation of [5], the 'convenient' gauge 
is equivalent to K = 0, G = 0, and j a (di) a = 0. A gauge is called good if an 
arbitrary perturbation p^ v can be transformed by a unique vector field to one that 
satisfies the gauge condition, p' = p^ v — V^H^ — V^S M . I call a gauge regular if 
p' v is everywhere regular whenever p M „ is everywhere regular. This is equivalent to 
saying that there must exist an everywhere regular vector field = (£ aj -C>A£) that 
implements the required transformation. Otherwise I call the gauge singular. From 



4 



appendix E of 5 , the explicit form of the gauge transformation is 

(8) 

(9) 
(10) 
(11) 
(12) 

(13) 
(14) 

The conditions G' — and j' a (d t ) a — j' a (d v ) a = j' v = are easily obtained by setting 
£ = r 2 G/2 and = j v — r 2 d v G/2, independent of any other requirements. On the 
other hand, setting K' = requires 
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Note that f(r — 2M) = 0, therefore, since K, G and j v were assumed to be arbi- 
trary regular functions, £ r must be singular at r = 2M, the future horizon, which is 
covered by our choice of advanced EF coordinates. Hence, the 'convenient' gauge is 
necessarily singular. Moreover, this singularity appears explicitly in the j' r , h' rr and 
h' vr components of the metric perturbation, though not in the h' vv one. 

Note that the above result would have been difficult, though not impossible, to 
obtain directly in static Schwarzschild coordinates. Since no point of the horizon is 
covered by (t, r) coordinates, even regular tensor fields may have components that 
diverge as powers of 1// as r — > 2AI, though such divergences will have a specific 
structure. For instance, this structure can be identified by transforming a regular ten- 
sor field from EF to Schwarzschild coordinates. Thus, the regularity of a tensor fields 
could be checked in static coordinates by examining the structure of the divergences 
of its components as r — > 2M, though with some effort. In fact, a variant of this 
technique was unsuccessfully, used in section 4.4 of [l0]. Unfortunately, the authors 
ultimately failed to notice the singular nature of their gauge choice. Also note that the 
same procedure can be used to check that the standard Regge- Wheeler gaug^ (j a = 0, 
G = 0) or the more recently proposed light-cone gauge of Preston & Poisson [6] are 
both regular. 

3.2 Zerilli function 



After imposing the 'convenient' gauge, the authors of 10 proceed to solve the con- 
straints among the remaining independent even-parity components of h^, isolate the 
single physical degree of freedom (the Zerilli function 4' Z ), an d express h^ u in terms of 
i\) Z . The Zerilli function is quantized as a normal scalar field. The expression for h^ u 
in terms of i\) Z is necessary to obtain the correct (lowest non-trivial order) coupling 



2 Referred to as Zerilli gauge in 101. 
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between i/) z and a Klein-Gordon field <f>. This coupling is obtained from the cubic term 
in the expansion of the standard, massless Einstein-Hilbert-Klein-Gordon Lagrangian. 

The derivation of the explicit expressions for in terms of ip z in 10 are fairly 
involved and specific to static Schwarzschild coordinates. However, given that the 
Zcrilli function is already known in gauge independent form, cf equation (4.23) of [5j, 
these expressions may be obtained algorithmically in any coordinate system, if the 
gauge-fixed equations of motion in that coordinate system are given. This opens the 
possibility of using computer algebra software to reduce the manual labor necessary 
to reproduce their derivation. 

For instance, the equations of motion for the even-parity metric perturbations in 
EF coordinates are given in appendix E of j5j. In 'convenient' gauge they reduce to 
the following. First, we write down the so-called gauge invariant combinations: 

J, -h ■ 2M f; „n 

"vv — "vv ~f ^2 JTl v ^/ 

h vr — h vr — v j r — 2~jn ("^-'v 

K = -—jr- (19) 

r 

In terms of them, the equations of motion for the metric perturbations, where for 
brevity we use A = 1(1 + 1) and \x— (£—!)(£ + 2), are 

= —d%K — -d r K — -d v h rr + —d r h rr + -d r h vr + —^-^—h r r, (20) 

TV T T 2r A 

= d v d r K + -d v K + - — ^-d r K — —d v h rr — -d r h vv (21) 

r 2 nvv 2r 2 r 2 2r 2 

r — M _ (r — M) f ~ 1 2 f f 

= -d 2 v K + —^d v K + V ^ U d r K + -d v h vv + -yd v h vr - J -8 r h vv (22) 

ur + 4M- _2/~ _p~ _/*/- 

2r 3 r 2 2r 2 

2 ~ r — M ~ 

Qyhyy dyhyp Qyl£ ~\~ Ijiyy \ ^ hyy ■ ^23^ 

= -0 v h vr + d r h vv - d v K - fd r K + 2{T ~ 2 M ' h vr + V—^ hrr, (24) 

~ r — M ~ 2 ~ 

= -d$hrr + 2d v d r h vr - d?h vv + 2d v d r K + fd%K -^—d v h rr + -d v K (25) 

2 - 2(r - M) - (r-M)f a - 2(r - M) - A - 

(Jyilyy (Jffb'D'p (Jyllyy ~\ 2 (JyJ\ 2 '"VT 



2 _ 

2r- hrr ' 



Xr 2 - 2uMr - AM 2 ~ 
h 



fh rr . (26) 



G 



The Zerilli function, from appendix E of j5j, is given by 



t 2r 



K + j(h rr - rd r K) 



(27) 



&. = (t-l)U + 2) + — =/i+— , ^>2. (28) 
r r 

The restriction on the spherical harmonic index stems from the fact that the I — 0, 1 
modes are not dynamical. The denominator A is clearly non-vanishing either in the 
exterior or interior black hole regions. Note that the above expression may differ by 
an ^-dependent multiplicative constant from ip z defined in |10| , which is immaterial 
for the purposes of this discussion. 

The algorithrrj^] for expressing p^ v in terms of \& consists of the following steps. 
Recall that each of the above equations is linear in the components of p^ and their 
partial derivatives. Also, let P n denote the set of all n-th partial derivatives of the p^ u 
components h vv , h vr , h rr , and j r . Note that we are not including ^> or its derivatives 
in P n . 

1. Initialization Let E be a list of expressions whose vanishing is equivalent to 



equations (201 to (27), e.g., the right hand sides of those equations. Further, 
divide this list into subsets E n , each containing no more than n partial deriva- 
tives acting on the components of p M „ (that is, variables from P Q up to P n 
only). Lastly, define E_\ be the subset of expressions that depend on ^ and its 
derivatives only; it starts out empty. 

2. Iteration Repeat for n = 1, 0, and — f , in that order: Apply d v and d r to each 
element of E n and collect the results in E' n+l . Using linear operations, eliminate 
the variables P n +i (being the highest order derivatives) from E n +\ U E' n+1 . 
Replace E n+ i by the eliminated expressions and add the remaining independent 
expressions to E n , which is possible since the remaining expressions will have no 
more than n derivatives acting on each component of p^ v . 

3. Termination Iterate step 2 until the number of independent expressions in Eq 
is the same as the number of variables in P . Optionally, keep iterating until 
E-\ is non-empty. 

4. Explicit Solution Set each expression in E to zero and solve the resulting 
linear equations for the variables in P . Each of the p^ v components will then 
be explicitly expressed in terms of \f r and its derivatives. 

If E-i is non-empty, then setting each of its elements to zero is equivalent to the 
explicit equation of motion for VP. This algorithm is not guaranteed to terminate for 
an arbitrary set of partial differential equations with constraints (though a generalized 
version of it is guaranteed to terminate under fairly general conditions). However, if 
it is known to terminate for a set of partial differential equations expressed in one 
coordinate system, then it will terminate for the same set of equations expressed in 
any other coordinate system. The results of section 2.1 of [To] essentially show that the 



algorithm terminates for equations ( 20 1 to ( 27 ) , when expressed in static Schwarzschild 
coordinates. 

3 This algorithm is inspired by the study of formal integrability of partial differential equations 
equations [s][9] . It is not too difficult to see how the simplified version presented here is equivalent 
to solving the constraint equations 'by hand.' 
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Applying this algorithm, we can obtain explicit expressions for the non-zero compo- 
nents of p^u in EF coordinates. The corresponding expressions in static Schwarzschild 
coordinates can be obtained by applying the same algorithm to the equations of motion 
explicitly given in appendix C of [5] or by applying the usual coordinate transformation 
rules of tensor calculus. The results agree with section 2.1 of [To] , 

In these explicit expressions, the singularity of the 'convenient' gauge is apparent 
from the presence of terms proportional to inverse powers of 1/ / in j r , h rr and h vr , 
which diverge as r — > 2M, presuming that is itself regular at the horizon. These 
divergences can be removed by the following explicit (singular) gauge transformation: 
£ = 0, = 0, and 

_ M (4McVg + AvE) 
er ~~2A ~f ■ (29) 

After this gauge transformation, the explicit expressions for the components of the 
metric perturbation become 



2/i*M(A + l) 2 (A + l)[12<9 r *M(A + 1) + /i*(A - 8) + 360„#M] 

3A 2 r 2 18Ar 
^>M 2 X + Ad v ^M 3 #M(A-8) /x* (A + 4) + 12 [i (d r $) M 

32, 



(30) 



+ r<9„* + d v ^ + d r ^ 
_2/j*M(A+l) //*(A + 4) - 12(d r *)M(A + l) - 18(9„*)M 
vr ~ AV^ 6Ar (31) 

#MA + 4 Af 2 3 (d„$) M + /x* 

2r 2 ~ ~~ 3r 

Ad^* + 2 (d 2 1<) (2M + r) - 2 (9„a r *) r + 4 (d r #) 
4 

6^Af 6(d r V)M 2(d v d r V)(2M + r) + 2(d v V)+ii(dry) 
" ~ A 2 r 2 Ar 2 ' (3 j 

3*M *A + 2 (d v V) (2M + r) + 2 (9 r *) r 
> = ^ 4 • (33) 

These expressions are manifestly regular at the horizon. Finally, the explicit equation 
of motion for \& is 



2d v d r V + fd^ + fdr* = ^ 



,2 



H + 2 6M\ 36M 2 / 2M\ 



/ r 4 \ r / 



*, (34) 



where the left hand side of the above equation is simply the d'Alambertian, on 
the (i, r)-plane. This is the well known equation of motion for the Zerilli function, 
derived both in 10 and (5l. 



4 Interaction and divergences 

Once both the scalar field and the metric perturbations are quantized, and their cubic 
coupling is explicitly derived, the calculation in [10] proceeds as follows (though this 
logic is only implicit in its technical details). The lowest order correction to the scalar 
2-point function is computed, which amounts to taking into account the single loop 
diagram shown in figure [I] This diagram is evaluated using an optical theorem-like 
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Figure 1: Lowest order correction to the scalar (straight lines) 2-point function from 
interaction with dynamical gravitons (wavy lines). The uj labels are short-hand for 
the mode indices and the sums may involve integrals over continuous labels. The loop 
diagram may be evaluated, using an optical theorem-like identity, by cutting it into 
two tree diagrams. 



identity, which is also schematically illustrated in figure. The correction to the scalar 
Hawking radiation spectrum is encoded in this correction to the 2-point function. 

On a Minkowski background, this kind of calculation is most conveniently per- 
formed in momentum space, where each leg of the interaction vertex appearing in 
figure [T] is labelled by an on-shell 4-momentum k and a tensor polarization index j 
(if any), while the vertex factor is proportional to the triple mode function overlap 
integral 

V k ,k>,k» = / d i x4>l{x)^ k ,{x)4>l, l {x) (35) 

= J d A x exp[-i(fc + k! + k")^] ~ 6(k + k' + k"). (36) 

It is apparent, since each mode function is everywhere regular, that the integrand 
defining the vertex factor is locally integrable, that is, its integral over any bounded 
region exists and is finite. It is also apparent that, nonetheless, Vk,k',h" is a distribu- 
tion, which follows from the global convergence properties of the above integral. Since 
these properties rely only on the local regularity of mode functions, they are expected 
to hold in curved space-times as well. 

On the black hole background, it is most convenient to perform diagrammatic 
calculations in angular momentum-frequency space. Each line of a diagram is then 
labelled by a frequency uj, a pair of spherical harmonics indices t, m, and a tensor 
polarization index j. Roughly speaking, the vertex factor is once again proportional 
to the triple mode function overlap integral 

Vvt^wu" = J tfi(r)^(t,r)^,(t,r)^,(*,r), (37) 

where w(r) is a 2-form on the (i,r)-plane, which takes into account the invariant 
volume measure and r-dependent coefficients that come from the expression for h^ u 
in terms of i/j z . Our discussion up to this point shows that, in a regular gauge, each 
of the w, 4>ui-i as weu as their products, should be locally integrable, including in the 
vicinity of the horizon. 



9 



More precisely, to take into account all kinds of vertices that couple the scalar 
field to metric perturbations, we must also consider derivative couplings. In that case, 



some of the <^w's in (37 1 will be acted upon by partial derivatives. If each ifi^e is 
regular, then any scalars made up of its derivatives will also be regular, therefore 
the local integrability argument is unmodified. Note, however, due to the presence 
of the black hole horizon, that the fail to be regular at the horizon. They are 
still bounded, but become highly oscillatory in the vicinity of the future horizon, 
with the oscillation phase diverging at the horizon itself (see section 10.2 of [2j, for 
instance). Hence, derivatives of (j>* t may become unbounded, though highly oscillatory, 
in a neighborhood of the horizon. Nonetheless, despite being unbounded, due to 
the oscillatory behavior of the integrands, their integrals should be evaluable in a 
distributional sense. Hence, even on a black hole background, and even with derivative 



couplings, the integrand in the triple mode function overlap (37) should be locally 
integrable, though perhaps only distributionally. 

4.1 Divergence in the triple mode function overlap 

The expressions (f° r even-parity modes) and (^ or odd-parity 

modes), introduced in equation (128) of [10] , are closely related to these kinds of 
triple mode function overl ap integrals. Despite the expectations expressed above, 



equations (160) to (162) of 10 show that the integrand defining K%' t ,. u „g, fails to be 
locally integrable. This failure of local integrability can be traced to the singularity 
of the 'convenient' gauge, which introduces terms of the form J dvol//" with n = 1, 2 



into the integral in (37), making it diverge in the vicinity of the horizon. In fact, as 
follows from the results of section [3j the only such locally non-integrable terms that 
contribute to the Ku'£'iv"£" are precisely the ones that can be removed by the explicit 



gauge transformation (29 1 and hence are pure gauge artifacts. On the other hand, the 
integrand defining V^^i^n would be locally integrable, as expected, in any regular 
gauge. 

Unfortunately, the authors of [lO] have mistakenly identified these divergent con- 
tributions to K^, as the dominant ones, have arbitrarily regulated them using 
a principal value prescription in the radial integral over an interval of size e about the 
horizon, and have dropped all other terms, including Hu'l'-u"l"- Since this e-regulator 
appears in the multiplicative coefficient of the spectrum correction in ([T]), the only 
conclusion to be drawn from the preceding discussion is that the size of the correction 
has not been correctly estimated and that what has been estimated is a pure gauge ar- 
tifact. To obtain a reliable estimate, the triple mode function overlap integrals would 
have to be analyzed anew, once they are rewritten in a regular gauge. 

4.2 Divergence in the summations over modes 



The logic outlined at the beginning of this section culminates in equation (144) of 10 
which expresses the correction to the Hawking radiation spectrum in terms of the triple 
mode function overlaps, K^,.^,,^, and H^,. u ,, e ,, the Bogolubov coefficients relating 
the <jyZtA modes, ot u ^ u i£i an d flue,u't'i an d some Clebsch-Gordan coefficients coming 



from the integration of products of spherical harmonics, Cp^,.^,^,. This expression 
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for the correction is schematically, as illustrated in figure [TJ 

(38) 



, „| 2 



where each compactly represents a combined sum-integral over all mode indices, 
also including spherical harmonic and polarization indices. Once V u ^',u" is estimated, 
the outer mode sums are seen to be divergent. In equations (163) and (164) of [Toj, this 
divergence is regulated by essentially introducing lower and upper frequency cutoffs, 



respectively 1/L and w cu t, cf also equation (82) in 10 



Leaving aside the fact that the estimates of the size of V^^' ^" cannot be completely 



trusted due to gauge artifacts, a divergence in (38) is to be expected, as in any generic 
1-loop perturbative calculation. The standard way to deal with such a divergence is to 
introduce a local counter-term in the original Lagrangian density. Since this divergence 
appears in a correction to the scalar self-energy, as illustrated in figure [ija) , such a 
counter-term would only renormalize the kinetic and mass parts of the scalar field 
Lagrangian density. Thus, renormalization would allow the regulator dependence of 
the final result to be removed, such that the coefficient in front of the correction in ([I]) 
would not depend on w cu t and L. 



5 Discussion 



The authors of 10 have tackled an interesting and challenging question. Unfortu- 
nately, their calculation suffers from a few problems. The so-called 'convenient' gauge 
chosen by the authors for the even-parity metric perturbations turns out to be singu- 
lar, unlike the standard Regge- Wheeler gauge. The singularities introduced by their 
choice of gauge result in spurious divergences, which mask all other contributions in 
the vertex factors characterizing the coupling of metric perturbations to the scalar 
field. Moreover, another divergence, corresponding to the expected 1-loop divergence 
of perturbative quantum field theory, is not removed via renormalization. Both kinds 
of divergences are regulated, introducing arbitrary parameters into the calculation. 
As clearly seen in ([lj, the final result for the correction to the spectrum of scalar 
Hawking radiation depends on regulators e, u; cut and L. Their presence, makes the 
given estimate for the size of the correction unreliable. 

Despite these problems, the authors have successfully addressed major, necessary 
parts of the calculation: (a) quantization of a scalar field and metric perturbations in 
a black hole background (via gauge fixing and explicit reduction to physical degrees 
of freedom), (b) explicit evaluation of the cubic scalar-graviton coupling (modulo 
gauge issues) , (c) estimation of the in-out mode Bogolubov coefficients via a geometric 
optics approximation and (d) an explicit expression for the spectrum correction in 
terms of Bogolubov coefficients, triple mode function overlaps and Clebsch-Gordan 
coefficients. It would suffice only minor modifications and a careful application of 
standard quantum field theoretic techniques to complete this calculation and obtain a 
definite, parameter-free estimate for the correction to the spectrum of scalar Hawking 
radiation. Moreover, the same techniques are readily applicable to the problem of 
perturbative back reaction of Hawking radiation on the quantum geometry of the 
black hole. 
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